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Correlated imaging through atmospheric turbulence is studied, and the analytical expressions describing 
turbulence effects on image resolution are derived. Compared with direct imaging, correlated imaging can 
reduce the influence of turbulence to a certain extent and reconstruct high-resolution images. The result is 
backed up by numerical simulations, in which turbulence-induced phase perturbations are simulated by random 
phase screens inserting propagation paths. © 2010 Optical Society of America 
OCIS codes: 270.0270, 010.1330, 110.0115 
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As correlated imaging develops well in recent years 
[1-5], more attention has been focused on how to ap- 
ply this technique to practical applications to over- 
come the limits in conventional optical systems. For 
an imaging system which must look through the atmo- 
sphere, turbulence-induced wavefront variations distort 
the point spread function (PSF) of the system from its 
ideal diffraction-limited shape, which leads to the the 
degradation of image resolution [6]. To mitigate tur- 
bulence effects, a number of methods, such as speckle 
imaging and adaptive optics techniques [6], have been 
proposed and applied in optical astronomy. Nonetheless, 
each of these techniques has its own set of performance 
limits, hardware and software requirements. New ap- 
proaches to the problem of reducing these effects are still 
of much interest. Here we investigate the performance of 
correlated imaging through atmospheric turbulence and 
find that the influence of turbulence can be weakened by 
the second-order intensity correlation. 

A schematic of correlated imaging through the atmo- 
sphere is depicted in Fig. 1. The beam splitter (BS) di- 
vides thermal light into two beams propagating through 
two distinct optical paths. One is test arm which includes 
an unknown object and a telescope setup consisting of a 
lens with focal length / and a detector Df. The object 
is located at a distance di from the source as well as d2 
to the telescope setup. The other is the reference arm 
where another telescope setup consisting of a lens and a 
detector Dr is placed at do = di -\- ^2 from the source. 
For remote sensing (i.e., di^d2 ^ /), the detector Df 
(or Dr) generally lies close to the back focal plane of the 
lens (i.e., ds ^ /)• The test arm is imbedded in the atmo- 
sphere, and turbulence- induced wavefront fluctuations in 
propagation paths di and d2 are represented by ^i and 
^2, respectively. While the reference arm is said to be 
a free-space propagation through the distance do by as- 
suming that there exists no turbulence. The assumption 
is based on the fact that the optical field in the reference 
arm is totally predictable if the field distribution of the 
source is well known [1,2]. 
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Fig. 1. Schematic of correlated imaging through atmo- 
spheric turbulence. 

In the test arm, the field Et{xt) in the detector Df can 
be given by 

Et{xt) = jj dxdiEs{x)hi{i,x)t{i)h2{xui), (1) 

where Es{x) corresponds to the source field, and t(^) de- 
notes the transmission function of the object. /ii(^,x), 
^2(^^,0 ^^^ ^h^ impulse response functions from the 
source to the object and from the object to the detector 
Dt^ respectively. 

Furthermore, according to the extended Huygens- 
Fresnel integral [7], hi{(^^ x) and /i2(^t, have the forms 



hi{£,,x) = 
h2{xt,C) 



^ ^^(^-0^+^1(^,0 
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where k = 27r/A is the wave number with A being the 
wavelength, and M = —ds/d2 is the magnification of 
the telescope setup. ^i(x, ^) and ^2(^,^) account for 
the random parts (due to atmospheric turbulence) of the 
complex phases of the fields in the propagation paths di 
and (^2, respectively. 



The field Erixr) in tiie detector Dr is connected to 
the source field Es {x) by the Fresnel diffraction integral 
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It's worth pointing out that the apertures of the lenses 
are regarded as large enough, and the diffraction limit of 
the lenses has been neglected here. 

Performing the intensity correlation measurement be- 
tween the test arm and the reference arm, we get 

G[Xt,Xr) = {It{Xt)Ir{Xr)) - {It{xt)) {Ir{Xr)) 

= Co I dxdx'dx"dx"'d(^d(,'{Es{x)El{x'")) 

x{E:{x')Es{x")){h,{(^,x)hl{i',x')) 

x{h2{xt.0K{xua)mt%a 
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(4) 



where cq is a constant {\^ did2d2,\M\)~^ ^ and 
It{xt)^ Ir{xr) represent the intensity distributions 
in Dt and Dr^ respectively. Here, we have supposed that 
the thermal field, and the two turbulent regions are 
statistically independent of each other. 

If the source is fully spatially incoherent and its inten- 
sity distribution is of the Gaussian type, the first-order 
correlation function of the source has the form 



{E,ix)E;ix'))=Ioe~ 



S{x — x')^ 



(5) 



where /q denotes the mean intensity at the center of the 
source, and r^ is the 1/e^ intensity radius. With the help 
of Eqs. (2a), (2b), and (5), Eq. (4) can be rewritten as 

G{xt,Xr) = /q I dxdxdr]dr]'didi't{i)e{i') 
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(6) 



The ensemble average of phase variations arising from 
turbulence can be approximated by [7] 



{e 



^i(a;,0+^*(^''^') 



_ -j^[{x-x'fHx-x'){e,-OH^-^'?] ,^, 

= e p^ , (1) 

where pi = {{)MbCl^^k'^di)-^^^ {i = 1,2) is the coher- 
ence length of a spherical wave propagating in the tur- 
bulent medium and Cn corresponds to the refractive- 
index structure constants describing the strength of at- 
mospheric turbulence in the propagation path di. It's 
worth emphasizing that we have adopted a quadratic 
approximation of the Rytov's phase structure function 
in Eq. (7) to obtain the analytical formula, and this ap- 
proximation has been used widely in literatures [4,7]. 



Substituting Eq. (7) to Eq. (6) and integrating over 
rj^T]' ^x^x' ^ we have 
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v/a/32(a + 2/3i) 
xe-^WT(^-^/^)%-f(^-*/^)',(8) 



where A = k/2di, B = k/2d2, a = r-'^/2, pi = p~'^ . 

By making Xr = Xt in Eq. (8), we carry out a spe- 
cial point-to-point intensity correlation [8] and obtain 
the PSF of the correlated imaging system 

h,{xr,0 = e-^^(«--/^)%-|^(«--/^)\ (9) 

For the sake of comparison, we also present the inten- 
sity distribution in Dt^ 



It{xt) 
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de|t(e)|'e-fe(«-^'/^)', (10) 



and the PSF of the test arm 



(11) 



From Eqs. (9) and (11), we can see that the full widths 
at half maximum (FWHM) of hg and ht both broaden 
with the increase of /3i (apart from the influence of 
the size of the source), which indicates that the resolu- 
tion, whether for correlated imaging or direct imaging, is 
degraded by atmospheric turbulence. Additionally, and 
most importantly, hg has a narrower FWHM compared 
to ht , which means that correlated imaging is helpful 
to reduce turbulent effects and achieve high-resolution 
images. 

In simulations, we consider correlated imaging through 
horizontal paths in the atmosphere, and thus C^ can be 
regarded as constant in the whole turbulent regions. The 
numerical model of light propagation in turbulence has 
been developed well [9,10]. The spatial power spectral 
density of the index of refraction fluctuations can be de- 
scribed by the Von Karman spectrum [9], 



^n{K, z) = 0m3Cl{z){K^ + L^2) 



-ll/6^-(K^o/27r)2 



(12) 



where K^ = K'^ -\- Ky -\- Ky ^ z is the propagation distance 
from the source, Lq and Iq represent the outer scale and 
inner scale of the turbulence, respectively. By using the 
spectrum in Eq. (12) to fllter a complex Gaussian pseu- 
dorandom fleld and inverse transforming the result, one 
obtains a two-dimensional phase screen which has the 
same statistics as the turbulence-induced phase varia- 
tions [9]. For long atmospheric paths, the multiple phase- 
screen model [10] has been used in simulations. The tur- 
bulent region with the propagation length di is broken 
into a number of layers with a thickness Az. Phase fluc- 
tuations in each layer are represented by a phase screen 
inserting in the middle of the layer. The effect of fleld 
propagation through these continuous layers can be cal- 
culated separately and then combined to characterize 
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Fig. 2. Simulated ( open circles) and theoretical (solid 
line) on-axis irradiance variance versus the propagation 
distance. The outer scale and inner scale of turbulence 
are Lq = 3 m and Iq = 1 cm, respectively. 



propagation through the entire turbulent region, pro- 
vided the index of refraction fluctuations for each layer 
are statistically independent [6]. 

First of all, to verify the computer programm, we in- 
vestigate the behavior of a Gaussian beam (waist ra- 
dius Wo = 7 cm and wavelength A = 2 fim) travel- 
ing through the atmosphere with a strong turbulence 
level (C^=10~^^m~^/^). The thickness of each layer is 
Az = 50 m. The on-axis normalized intensity variance, 
defined as aj = (/^)/(/)^ — 1 [9], is plotted as a function 
of the propagation distance in Fig. 2. The good coinci- 
dence between the simulated data (open circles) and the 
theoretical result (solid line) predicted by [11] proves the 
validity of the programm. 

After the validation, we apply the programm to simu- 
late the correlated imaging system shown in Fig. 1. The 
thermal source (A = 0.532 /im and diameter D = 2re = 5 
cm) was described by a grid of 512 x 512 with a sample 
spacing Ax = Ay = 5 mm. The distances were set as 
di = d2 = 10 km and the focal length / = 1 m. The 
turbulence regions in the paths di{i = 1,2) were divided 
into 20 layers with a thickness Az = 500 m, respec- 
tively. The turbulent parameters were assumed constant 
at the outer scale Lq = 100 m and the inner scale Iq = 5 
mm. By averaging over 10^ samples, simulated results 
[see Fig. (3)] clearly show the image resolution decrease 
with the increase of the turbulent strength, which ac- 
cords with the analytical calculation from Eq. (8). 
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Fig. 3. The reconstructed images (from left to right) via 
the correlation in the atmosphere with turbulent levels 
C^ = lO-i^m-2/3 2.5 X lO-i^m-2/3 5 x lO-^^m-2/3 
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To compare direct imaging and correlated imaging, a 
simple double slit (slit width 10 cm and center-to-center 
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Fig. 4. The acquired images of the double slit in the 
atmosphere with turbulent level C^ = 10~^^m~^/^ . (a) 
was obtained by the test arm, and (b) was extracted from 
the correlation. The normalized horizontal sections of the 
images are plotted in (c), where open circles correspond 
to the simulated data and solid lines show the theoretical 
predictions from Eqs. (8) and (10), respectively. 



separation 20 cm) was used. After statistics over 10^ 
samples, we obtained a blurred image detected by the 
test arm directly [see Fig. 4(a)] and a clear image re- 
constructed through the correlation [see Fig. 4(b)]. This 
confirms the analytical result that ghost imaging could 
reduce turbulent effects and improve resolution. 

In summary, by taking advantage of the extended 
Huygens-Fresnel integral, we have presented the theoret- 
ical expressions that describes how atmospheric turbu- 
lence corrupts the image resolution. Meanwhile, the an- 
alytical calculations and the numerical simulations have 
demonstrated that correlated imaging can provide imag- 
ing performance superior to direct imaging through the 
atmosphere. As an unique image-formed method, cor- 
related imaging can be effectively combined with con- 
ventional phase compensating techniques (e.g., adaptive 
optics) to further eliminate turbulent effects. 
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